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Abstract 


We  consider  the  problem  of  interpolating  a function  on  the  unit 

square  from  values  of  the  function  given  on  a set  (called  the  "^^sign^T^ 

of  N discrete  points.  -M^ot>t»4t>^ome  new  results^on  better  ways  to 

distribute  these  N points  so  as  to  minimize  the  maximum  interpolation 

error.  Gordon  [8,  Approximation  With  Special  Emphasis  on  Spline 

Functions,  I.  Schoenberg,  ed.  1969],  showed  that  an  improvement  on  the 

tensor  product  design  {(x,y)}  = (^  , ^ , i , j = l,2,.,.,n},  of  N = n 

interpolation  points  could  be  made  by  using  a design  that  is  the  union 

of  two  tensor  product  designs:  {(x,y)}  “ i “ 1,2,. ..,n, 

j = 1,2,. ..,n  } u {-y  , ^ , i = 1,2,. ..,n  , j = 1,2,. ..,n}.  Here 
o 2 

N = 2 n - n since  the  intersection  of  these  tensor  product  designs 
has  n points.  Gordon  showed  that  where  interpolation  on  the  original 
tensor  product  design  has  a convergence  rate  0(N  for  some  p with 
respect  to  functions  with  specified  smoothness  properties,  the  new 
design  will  have  an  associated  convergence  rate  of  We 

construct  the  tth  design  which  is  the  union  of  t+1  tensor  product  designs: 
{(x.y)}  = U ^ . i = l,2,...,n'^,  j = l,2,...,n  . Here 

, 0+2  )l+l 

the  tth  design  has  M = ()i+l)n  - tn  distinct  points,  and  we  show 

that  the  associated  convergence  rate  is  Some  exact 

theorems  to  this  effect  are  proved  for  convergence  of  minimal  norm 
interpolation  in  tensor  product  reproducing  kernel  Hilbert  spaces.  An 
application  to  the  reduction  of  X-ray  dosage  in  the  reconstruction  of 
human  internal  structure  from  X-ray  projections  is  suggested, 

Applications  to  the  desion  of  certain  (other)  indirect  sensing 
experiments  are  also  noted. 


1.  Introduction 


We  consider  the  problem  of  interpolating  a surface  on  the  unit 
square  from  values  of  the  surface  given  on  a set  (to  be  called  the 
“design")  of  N discrete  points.  We  obtain  some  new  results  on  better 
ways  to  distribute  these  N points  so  as  to  minimize  the  interpolation 
error. 

Let  f(x,y)  be  a function  defined  on  [0,l]x[0,l]  with  the  property 

that 

— _ f(x,y)  exists  and  is  continuous  (1.1 

ax^ay^ 

for  0 _<  p <_  2m  and  0 £ q £ 2m.  Convergence  properties  of  piecewise 
polynomial,  including  bivariate  spline, interpolation  are  known  for 
such  functions  when  interpolation  is  carried  out  over  a tensor  product 
design.  By  a tensor  product  design  is  meant  a set  of  N = N1N2 
points  in  the  unit  square  of  the  form  {(x^,yj)  i = 1,2,...,N^, 

j = 1,2 N2}  where  0 £ x.|  < ><2  ^ ^ ® - ^1  ^N2  - ^ ' 

See  deBoor  and  Fix  [6]  Lyche  and  Schumaker  [14J.  Specializing  to 

Ni  = N2  = n,  N = n , in  general  one  has,  with  appropriate  interpolants  f, 

sup  lf{x,y)-f(x,y)|  = 0(-^)  = O(^) 
x,y  n N 

as  n -►  ®.  This  is  the  same  rate  as  is  achieveable  in  univariate 

Interpolation  - if  g(z),  z e [0,1]  has  2m  continuous  derivatives,  then, 

for  example,  local  Lagrange  (polynomial)  interpolation  of  degree  2m-l 

through  2m  adjacent  points  ^ ....  gives  a remainder  term 

0(-^),  see  [12,  p.l90],  and  the  same  rates  are  obtainable  for  various 
n 

types  of  spline  interpolation  (see,  e.g.  [22]). 
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It  will  :o  convenient  here  to  report  our  results  for  functions  in 
a reproducing  kernel  Hilbert  space  (r.k.h.s.)  of  functions  which 
satisfy  (1.1)  for  p,q  £ 2m-l  and  the  slightly  weaker  condition 


1 1 


,4m 


/ ^ ^~2m~2i  ^ "• 

Convergence  results  for  certain  types  of  bivariate  spline 
interpolation  on  the  tensor  product  design  q ^ ^n  * ^^i  j=l  * 

N = n^,  for  functions  in  an  r.k.h.s.,  satisfying  (1.2)  have  been 
obtained  by  Mansfield  [17,  Eq.  (4.2)],  who  gives 

sup  |f(x,y)-f(x,y)|  = • 

x,ye[0,l]  n ' N ' 


(1.2) 


and  it  appears  that  this  is  the  best  possible  rate  under  the  given 
conditions,  for  the  tensor  product  design  q. 

For  many  years  we  have  been  intrigued  by  W.  J.  Gordon's  results 
[8]  with  the  so-called  blending  function  splines.  Gordon  used  designs 
of  the  form 


N.l 


where 


jJ  * k _ ^ C n -1  n 


(1.3) 


1 2 

The  two  tensor  product  designs  T ’ and 


2 1 3 

’ each  have  n points  and 


their  intersection,  T^*^  0 Tp'^’  ^n*^  " ^n  ’ ^^c=l  ,n=l  ’ 

2 3 2 

n points,  so  that  their  union  has  N = 2 n - n distinct 

* 12  2 1 

points.  See  Figure  1.1,  which  depicts  the  designs  Tp’  , T^’  and  Tj^ 
for  n = 5, (thus  N=225).  Figure  1.2,  taken  from  Gordon,  schematically 
describes  the  interpolation  of  f over  part  of  T^^j. 

/ 


Gordon  showed 


-3 


^225,1 


Figure  1 .1 
The  designs  and 


n=5,  N=225 


-4- 


theoretically  and  by  numerical  example  with  bicubic  spline  interpolation, 
that  interpolation  to  appropriately  smooth  functions  on  the  N points 
distributed  according  to  the  design  -j  gave  a better  approximation 
than  interpolation  on  the  N points  distributed  according  to  the  tensor 
product  design  q. 

Letting  Qf  and  be  points  in  and  Gordon 

showed,  loosely  speaking,  that  if 


sup  lf(x,y)-{P^ 

x,y  nP 

for  some  p,  then 

sup  lf(x,y)-(Pj^^^f)(x,y)|  = 0(-^)  = 0(^g-p-^^)  . (1.5) 

(Here  we  are  using  the  fact  that  n = , i = 0,1.) 

In  practice,  if  N is  large,  the  difference  between  (1.4)  and  (1.5) 
can  be  important.  Recently  Ylvisaker  [28]  gave  a similar  theorem 
comparing  convergence  rates  with  T|^  q and  i the  context  of 
minimal  norm  interpolation  in  an  r.k.h.s.  He  observes  that 
Improvements  over  -j  should  exist,  but  apparently  are  not  available. 


Spurred  on  by  Ylvisaker's  existence  conjecture,  we  began  a search 
for  designs  that  improve  upon  the  rate  0( '2'p/ j) • search  was 

successful,  and  in  this  paper  we  provide  a sequence  of  designs 
Tu  .•  I - 0,1,...  with  the  following  property.  When  the  tensor  product 
design  q results  in  a convergence  rate  of  0(-^)  = 0( -py^l  t the  design 
Tfj  ^ with  i fixed,  has  N = (t+l)n  -tn*'  points,  and  results  in  the 
convergence  rate 
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as  n Considerable  improvement  in  some  practical  problems  appears 

to  be  possible  as  a result. 

Later  in  this  paper  we  suggest  how  this  result  might  be 
used  to  reduce  the  X-ray  dosage  in  computerized  image  reconstruction 
while  maintaining  picture  quality.  Other  potential  applications  which 
come  to  mind,  are  patterns  of  weather  balloon  releases  or  other 
weather  measurements,  and  patterms  of  seismic  soundings,  and  ocean 
depth  measurements. 

The  designs  we  have  found  take  the  following  form: 

j k 

Let  be  the  tensor  product  design  defined  in 

(1.3)  and  let 

T.  , = u I = 0,1,2,...  (1.6) 


The  1=0  and  i = 1 designs  have  already  been  described, 
shows  tJ’^,  T^’^,  T^’^  and  T^  2 ^ = 189. 

It  is  easy  to  see  that 


•rjjk+r  p.  •|.j+s,k  _ 

^n  ''  'n  ■ 'n  ’ 


Figure  1 .3 


for  r,s  non-negative  integers,  and  it  can  be  seen,  with  reference  to 
the  Venn  diagram  of  Figure  1.4a,  that 
1+1  1 


T - V Tj»l+2“j  y -pjjl+l'j 


(1.7) 
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Figure  1.3 

The  designs  T^’^,  T^’^,  and  for  n=3,  Ji*2,  N=189 


X 


Si 
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I 


i k 

In  Figure  1.4,  ^ is  outlined  boldly,  and  each  set  Tj^’  is 

represented  by  the  rectangle  with  upper  right  hand  corner  in  position 
j,k.  The  Venn  diagram  for  -j  is  given  in  Figure  1.4b.  Since 

l-j.t+2-j  (contains  points  and  contains  n*"^^  points, 

n+p  o+l 

T|^  contains  N = (Ji+l)n  -tn  points. 

The  exact  convergence  theorems  we  prove  are  for  minimal  norm 
interpolation  in  a tensor  product  reproducing  kernel  Hilbert  space 
of  functions  on  [0,l]x[0,l],  to  be  called  Pi^operties  of 

tensor  product  r.k.h.s.  :hat  we  need  will  be  reviewed  in  Section  3. 


Let  be  the  (Sobolev)  space  of  functions  on  [0,1]: 

= {g,:  q,g; . . . ^ abs.  cont.,  g^'^^e  L2[0,1]}.  Then  we  let 

^K(m)  ~ ® or,  slightly  more  generally,  ® ^q* 

where  Hq  is  an  r.k.h.s.  of  functions  on  [0,1]  possessing  property  m to 


be  defined.  Functions  in 


will  satisfy 


gP+q 

— - — - f(x,y)  continuous,  P.  q 1 m-1 
3x^3y^ 

1 1 „2m  « 

f f mHii  ^ “• 

0 0 3x'"3y'" 

If  (for  example)  is  endowed  with  the  norm 

llgll^mi  = V (q^^’^(o))^  + / (g^'^^u))^  du  , (i.s) 

j=0  0 

then  minimal  norm  interpolation  in  is  polynomial  spline  inter- 

polation of  degree  2m-l  (see  [14J)  and  minimal  norm  interpolation  in 
0 can  be  shown  to  be  a form  of  bi-(2m-1)-ic  polynomial 
spline  interpolation.  We  argue  that  our  convergence  results  must 


hold  more  generally  for  bivariate  spline  interpolation  (that  is, 
not  necessarily  minimum-norm,  usually  a semi-norm  is  minimized)  but 
the  proofs  are  omitted.  Our  main  result  is  as  follows. 


Theorem 


Let  0 (more  generally,  let  ^ ® 

where  Q possess  "property  m"),  let  Pj^  ^ be  the  minimal  norm 

interpolation  operator  to  the  points  in  the  design  Tj^  ^ given  by  (1.6), 

0+1 

hence  N = (ii+l)n  -in  . Let  ll’ll|^  be  the  norm  in  and 

let  f e H|^.  Then,  for  each  fixed  i - 0,1 


' y€[0  - (n<^+l^Ml/2)  1 1^1 

(1.9) 


7S(M1/2)) 

= 0[()i+2)(-y-)^  ^ ] , 


(1.10) 


as  n “ 


where  c and  d are  two  constants  depending  on  Q.  If,  (loosely) 
furthermore  f e M|^j2ni)  satisfies  some  boundary  conditions,  then 
as  n 

2 £+1  |S(2ni-(l/2)) 

sup  |f(x,y)-(P.  /)(x,y)|  = 0( (ii+2)^(^)^^2  ).(1.11) 

x.ye[0,l]  ^ 

The  £ = 0 results  correspond  to  known  results  for  interpolation 
on  tensor  product  grids,  in  particular  (1.11)  corresponds  closely  to 
Mansfield  [17],  equation  (4.2). 

To  get  a feel  for  numbers,  consider  a typical  tensor  product 
grid  of  N = nxn  = 200x200  = 40,000  such  as  might  occur  in  reconstruction 
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of  the  internal  structure  of  the  human  body  from  multiple  X-ray 
projections.  X-ray  dosage  is  proportional  to  N.  (See  Section  4 
for  further  details).  In  the  application  described  in  Section  4 
it  is  argued  that  m = 1 is  appropriate,  (1.9)  holds  and  c = d. 

r 

O 

Letting  c = d,  e = c ||f|||^  in  (1.9),  the  error  bound  is 
"li+Vy/?  ^ ■ 0,1,...,  we  obtain  the  error  bound  for 

the  tensor  product  grid  q as 


200^/^ 


.1420  . 


£ 

Consider  alternatively  the  design  with  n = 10,  «,  = 2, 
N = (t+1  )n*'^^-An^^^  = 28,000.  The  error  bound  becomes 


e_ 

1,000 


TTI  " 


.1360  . 


Thus,  one  can  maintain  about  the  same  error  bound  with  a 30% 
reduction  in  N.  One  of  the  Gordon  designs  T^  i is  better  than  the 
tensor  product  design  for  example  = 1 , n = 27  gives  N = 39,366, 
with  a reduction  of  the  error  bound  to  .1140. 


In  applications  where  the  function  to  be  reconstructed  has  a 
greater  degree  of  smoothness,  (i.e.  m > 1)  the  relative  improvement 
can  be  very  much  greater. 

In  Section  2,  we  discuss  interpolation  operators  of  a fairly 
general  class,  and  obtain  formulas  for  the  remainder  operator  I " Pjj 
where  ^ is  any  interpolation  operator  of  a general  class,  not 
necessarily  minimum  norjn.  In  Section  3,  we  review  certain  properties 
of  tensor  product  r.k.h.s.  and  prove  (1.9),  (1.10)  and  (1.11).  In  Section  4, 
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we  suggest  the  possible  application  of  the  results  to  reducing  the 
X-ray  dosage  in  computerized  image  reconstruction.  In  Section  5 
we  mention  another  potential  application  and  some  related  unsolved 
problems,  and  in  Section  6,  we  indicate  how  certain  design  problems  for 
interpolation  are  equivalent  to  desinn  problems  for  indirect  sensing 
experiments. 

2.  Interpolation  Operators.  Remainder  Formulae 

In  this  section  we  develop  properties  of  a class  of  interpolation 
operators.  The  results  extend  those  of  Gordon  [8,9]. 

M M 

Let  be  any  M real  valued  functions  on  [0,1]  with  the 

property 


^ o’  ] ! j , i,j  = 1,2,...,M 


(2.1) 


For  f a real-valued  function  on  [0,l]x[0,l],  define  the  operators  P , 


and  Py  by 


Ni  '^l  N,  . 

(Px  f)(x.y)  = I <}>.•  (x)f(jr-,y) 

N ^2  N 

(p/f)(x,y)  = I <{i/(y)f(x,ij-) 

^ j=l  J 2 


(2.2) 

(2.3) 


and 


^1  ^2  N, 


(P  ' ^f)(x,y)  = I I ^ >(kU  Uy)f(^,^) 

i=l  j=l  ’ J *^1  *^2 


(2.4) 


It  is  obvious  that 

N,  N.  N„  N,  N,,N5 

Px‘  Py  = Py  Px  - P • (2.5) 

N,  ,N2 

and  that  P f interpolates  to  f for  all  (x,y)  on  the  tensor  product 
. Ni  N, 


i 1 ”1  ’2 

design 


The  reproducing  kernel  (r.k.)  R(o,t).  o,  i c T.  for  the  space  H,  is 
then  defined  by 


± 


4i 
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y(t)  • g(tKc{t],  t ■ Mhere  c(t)  Is  a random  variable  with 


w o 
o a* 


^ «»■ 
a ^ 
rsj 


o c 
o 

a.  4-» 

4Q  u 


I 

1/  hx. 
U 


U 

351 


N JZ 

w 

a>  s. 
.*-> 

t>  • 
IQ  V. 
k-  at 
IQ  Lr 
JZ  Zi 

u sz 


3 < 
JD  cn 
oi  at  • 


»—  o iq' 

ro  -M  £■ 

S;^  . 

COO 

S Sr: 

V.  < 
• £!.►-■ 
IQ  a</> 


o o CSJ 

^ IQ  rx. 

at  a 

W1  ^ 


at  tQ 
c J=  u o 
at 

o > c 


O • • 

5 I 

<0  O'  »0  *3 

■o  at  -o  «- 

■a  'p-  tQ 

>>  r-  a 

tA  iQ  ^ 

o > Q 

0 at  I 

c w 4/t  «/) 

V.  l/>  o 

01  o p ••- 

c o C 4J 

•p-  W kt) 

£ at  -p-- 

.c  o ■«  rx 

o at  IQ  r*x 

0 rsi  4-t  Ot 

1 g'--'" 

•p-  IQ  C U 

• 4-t  k.  o at 
iQ  « at  kA  a 

2 E c - o 
^ at  tj  *-» 

19  CTV^  O 


2> 


o u a 
o to  a 
T9  at 
o c • 

L.  -f-  «r 


u 
«Q  at 
u > 


8 ^ at 
H- 

to  • 
*A  • -kt  ^x 

at  <Q  00  ^x 

Otp—  O' 

gi^"- 

C at  •> 


o u ^ 
kk-  »*-  c u 
o at  o 

C 'p*  >■ 

coo 


o > o ^ 

O C E 

at  o k*-  < 
o o o 

IQ  • 

■o  OS  -M  o 
at  c f— 
N c at  to 
•P.  o E *k- 
O 

at  T9  o 3 

■M  at  to  03 
3 «/>  a 


lA  >t  • 4-t 


•<-■  E to 


T3  at  o p- 

•p-  rx  a X k- 

^ CD  O'  <A  o at 

tQ  P-  U M 

o -c  >1  a t-t 

•p-  4-t  p ^ a 3 

4^  rsj  tA  < CD 


C LTt 

at  cvj 

»—  O' 


<4*  O 
p O 

s_  •— 


a IQ  <y  ^ 

E a o IQ  >.  I 


k.  at 
p at  a 
T?  o 


03  tA  • JX 
C Ax  I- 

• O O'  o 


• o a 
-o  u at  k*. 
c aos  o 


O 4-t 

to 

E C 
c ‘P-  IQ 

tQ  X "P- 

• £.  IQ 

a c 
a IQ 

p to  o 


— E — 
o p 
a C (A 
o V*  o> 


C «A  TJ 
•p-  at  at 

•p-  ct 

C 4-t 

o t.  o 
at  u. 
p aa 

S at 
. a u 
•-*  c 

»-  01 
P to  k. 

8*  ^ »2^ 
at  c c 
a o 

CEO 


i.  oo 
O Ax 
•4-  O' 


k-  o 

k.  I- 


s 


CM 

* ^ 


fe  S 

(A  lO 
C Ot 

at  kA 

"5 


ig 

to  at 
^ a 

3 

a^ 


ar  «A 


! . •SS  . 
it  tit 

O'  o O' 
m ^ Lt  «A  ^ 
> C 

. • 4^  o U 

cn  c u 01 
L>  at  kA  a 

9)  p -p.  -p-  e 

3 ^ u 3 at 

< I—  -p*  > 


o c • 

C3  -P-  kO 

^ lA  V 

T)  «A  O C 

c c a o ^ 

to  o 4-*  u 

at  «A  44 

. g - 1. 

a u ^ P 

c at  a 

• 3 a H-  at 

e k4>  44  o o: 


o • a ^ 

<M  to 

PAX.  P e 

«-  O'  P-  < 

i.4-  i. 


iTI 

o ko  p as 

s & 


P 44  O 
‘P*  ^ to  P 
k kA  dO 


a CD 

rx 

c o> 
o «- 

44  P 

g- 


S£1 

iQ  O 


a x| 

kA»  g! 


O kA 

■o  c 


«A  m 
C CD 

a CM 


p S» 


c.  g 

p a 

oe  e 
pS 


c u 
® a * 
a to  T3 


o S u 
a kA  ^ 
a c c 
to  o at 

.p-  o 

at  4t  f 

a Q u 

44  E (/) 

■o  x . 

c p PQ  . 

IQ  k.  .CO 
a "3  rx. 
kA  a CM 

at  < p I 

0 3 4A  CD 

•P-  C CM 
44  C o CM 

iQ  44  a 
p-  * u a 
kA  c 

at  c 3 p 
> o a O' 

•fp a 
♦J  44  at  O' 
3 O C P- 

a c .pp 

3 P—  p 
w k4.  a kA 
44  kD  kA 

kA  at  at 

^ 44  c t- 
o iQ  o a 
p 1.  lA  o 

. A»  >pp  •<» 

•^  > kA  E 
• IQ  Ot 
:»  44  £ xi 
p—  a <Q 


Si 


•P-  c 

k.  o 

to  -P- 


44  to  O' 

kA  4-t  ^ 

at  kD 


at  "o  3 


4-»  44 

kA  U 

k.  •»  O 
o :x  a 

<4-  01 

«t-  a 

E 

I >,  iQ 


at  >t  o 
•—  44  a 
a Ot 


a at  «— 

> to 

at  o 
a c -p- 


ot  Wt  44 

c p e 

•p-  C p 


E c a 

kAP.^ 

at  kA 

C U kA 

z e.ie 

a u 44 

kA  at  kA 


o ft- 

>»S 


■ 44  O ^ ^ ^ 


•p-  > at 
c — r- 
o c 

C 3 p 


(A  at 
< at  E 

O *4 
p to  k. 
. •*-  IQ 

CO  k.  a 

.35 


a c a 
e Ok  f-  at 
•PI  c kA  a 
c 

p ^ o p 
• k-  kj  m 
CO  to  kA  o 


kA  i.  « 
c to  rx. 

at  01  o 


V (A 
PC  k 

• to  a 


t.  -P-  u% 

at  «A  ^ 


t-  >>fD 

44  P—  kA 

<J  -O  I 


U • 
>'rx. 
u * 
o rx. 
44  so 
IQ  • 


gx 


CT  at  I— 

c *4  rx. 
•*-  C O' 
T3  p— 


kA  acD 
iQ  aCD 
a to 


at  3 44 
44  E to 

kA  U 9r 

k.  o 

•o  a . 

iQ  c eii 
:z  o aJ 


e I 

cS| 

«*>  a 
«t  a 
a 

s;.-- 

0 a a r^ 

1 . .82 

at  •—  u 

44  fkk  h-  • 

52  c-- 

...  . °« 

O » M CD 

« - 2 r 

OOO  3 l>| 


^ C3 
IQ 

kpt  O 


‘CD  • *- 


fc.  »- 

at  t«| 

Si' 


44  3 


-i 


33 


c:  8 k^  < c ^ ^ 


C P- 

o :r 


p ^ o S 
o 8 


p^  k I 


S 


1 


^ s 


^ CM 


ID  « U) 

u a4  t— I 


r*  CM 


iCCMMiTT  CbAUirtCATiow  tmi*  ^agc  fwrnm  dm*  r« 


